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1^ ■ Abstract 

The expected areas of the Wiener sausages swept by a disc attached to the two- 
dimensional Brownian Bridge joining the origin to a point x over a time interval 
[0, t\ are computed. It is proved that the leading term of the expectation is given by 
^ ■ Ramanujan's function if |x| = 0{y/t). The second term is also given explicitly when 

CN . |x| = o{^/t). The corresponding result for unconditioned process is also obtained. 

d 



1 Introduction 

Let Bt be the standard two-dimensional Brownian motion started at the origin defined 
on a probabihty space (fi, J^, P). Fixing r > let 5'^^'''' be the Wiener sausage of radius r 
and length t, namely it is the region swept by the disc of radius r attached to at its 



. 5t ■ center as s runs from to t: 

5f ^ = {z G : - z| < r for some s G [0, t]}. 

In this paper we compute the expectation of the area of Si~\ which we denote by 
Area(S'j^'*), for Brownian motion conditioned to be at a prescribed point x G at time 
t as well as for free (unconditioned) Brownian motion. Define A^(A), called Ramanujan's 
function ([2], [21]) or integral ([7], page 219), by 

JO (Ig UY -|- TT^ U 

Put K = 26"^"^ where 7 = — e~^\gudu (Euler's constant). Let E designate the 
expectation with respect to P. In this paper we prove the following two theorems. 
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Theorem 1.1 



|E[Area(^f))] = 2vriv(^) + ^^^^^(1 + o(l)) as t ^ oo, 

— + ^ + 0(Vi) as t^O. 



t 2 

For X G we write x^ for the square of Euchdian length |x|. 
Theorem 1.2 For each M > 1, uniformly for |x| < My/t, as t —> oo 

' nt\ TTX^ r, / t 



X2 V 1 



0(11 



+ 0{l] 



where the left-hand side is the conditional expectation conditioned on Bt = 

The proof of Theorem 11.11 is performed by Laplace inversion with rather simple com- 
putations. The leading term in the formula of Theorem 11.21 is derived in a similar way to 
that in Theorem II. H while the identification of the second term of it is made by a delicate 
analysis unless |x| remains in a bounded set. 

Remark 1. The function N[t) admits the following asymptotic expansion in powers 
of 1/ Igt as t — 7- oo: 

N(f\ 1 "< 7' - C(2) , 3iC(2) + (-1? - 2C(3) 

where C(^) = Z^^i '^"^^ While it requires a tedious computation to derive the corre- 
sponding expansion of N{nt/r'^) directly from this one, actually there is a simple way to 
transform the expansion. For each a > the expansion for N{at) is obtained by simply 
replacing —7 by —7 — Iga in ( II. ip . so that the expansion up to the third order term 
becomes 

N(r.,\ ^ -7-lga , (^ + lga)^-C(2) ^ 

N{at) — + ——— + — — + ••• 1.2 

Igt (Igt)^ (Igt)'^ 

Similarly the asymptotic expansion of | Jq N(as)ds is obtained by replacing —7 by 1 — 
7 — Iga and — C(^) by 1 — ({k) in (11. ip . so that 

1 f .r. . , 1 ^ l-7-lg« ^ (l-7-lga)2 + l-C(2) 

-/ N{as)dsr^ — + — — ^ + — — ^ + ■■■ 1.3 

t Jo Igt [Igty (lgt)3 

(See Section 5 for these and related matters.) It also is noted that combining Theorems 
O and O yields 

E[Area(Sf ^)] - ^[Area(5f ^) \ Bt = 0] = 27i t\- asN'{as)\ds + 0(1) 

Jo L J 

as t —7- 00, where a = /t/r^. The asymptotic expansion of this difference is obtained from 
that of E[Area(Sf ^)]/lgt by a very simple rule (see Lemma [5.21 in Section 5). 

Remark 2. It may be reasonable to compare the expected increment of the sausage 
at time t with 27rr£'[|i?(|] as t J, 0. The former one is asymptotic to 2r\/27it according to 



the second half of Theorem ll.il while the latter equals 27rry 7rt/2, so that the ratio of the 
latter to the former equals 7r/2. 

Asymptotic behavior of the area Area(5'(^'^^) (or volume in higher dimensions) as t — )■ cxd 
has long been studied from various points of view. It is a typical functional of Brownian 
paths that is non-Markovian and the standard limit theorems for it (the law of large 
numbers, the central limit theorems or the large deviations) have been of continued interest 
([30]. [16], [To], etc.) . The expectation of it is the total heat emitted in the time interval 
[0, t] from the disc which is kept at the unit temperature. The sausage for conditioned 
process (with x = 0) naturally arises in the study of the asymptotic estimate of the 
trace of the heat kernel on the plane with randomly scattered cooling discs kept at zero 
temperature and has been effectively used (cf. [13j, [5j). 

For free Brownian motion, some asymptotic expansions in negative powers of Igt of 
the expectation divided by t for the sausage swept by an arbitrary (non-polar) compact 
subset of are obtained by Spitzer [22] up to magnitude o(l/(lgt)^) and by Le Gall 
[Is] to any order (see also [12]). M. van den Berg and E. Bolthausen [3] compute the 
expectation for Brownian motion conditioned on returning to the origin at time t for the 
disc case up to the error term of magnitude 0{t^J\g\gt)/{\gt)'^ and conjecture that their 
formula for the conditional expansion would be valid for the general compact set, K say, 
of positive capacity if the radius of the disc is replaced by the logarithmic capacity of 
K (in the usual normalization [T]), the situation already observed for the free Brownian 
case. 

The conjecture stated in [3J has been verified in a very recent paper [L9l by I. McGillivray: 
in fact he obtains the asymptotic expansion for any compact set K of positive capacity 
and computes the explicit forms of the first three coefficients of the expansion in terms 
of the logarithmic capacity of K, the coefficients agreeing with those conjectured in [3]. 
It is warned that these authors define the sausage for the Brownian motion Bt = i?2t 
instead of Bt so that to translate our results to their case one must replace t by 2t in our 
formulae. The coefficients of our formulae of course agree with those obtained previously, 
whether it is a free or conditioned Brownian motion which forms the sausage, as readily 
ascertained by substituting a = k/v"^ in (11. 2p and (II. 3p and noting that the logarithmic 
capacity of the disc of radius r equals (according to the normalization of logarithmic 
capacity in them). This suggests that our formulae in Theorems 11.11 and 11.21 would be 
extended to the sausage swept by any non-polar compact set K in place of the disc if r is 
replaced by the logarithmic capacity of K. 

The higher dimensional case of free Brownian motion is treated by Spitzer [22], [8] 
and Le Gall [15] for a compact set and by [9] for a ball. The pinned cases are dealt with 
by Uhlenbeck and Beth [28] and McGillivray [18]. A more detailed account of the results 
for the expected volume (for dimensions > 2) obtained up to 1997 can be found in [3]. 

The function N{t) expressing the leading terms in our theorems has already appeared 
in an analogous manner for the corresponding problem concerning the range of random 
walks (cf. [21]) (in which N'{X) is denoted by W{\)). The expression by means of N(t) 
makes possible a subtle comparison between the expected area of the Wiener sausage and 
the corresponding quantity to a random walk of mean zero: we can associate a certain 
natural 'radius', say r^,, of a lattice point of with the random walk ([25], Remark 6), 
and according to Corollary 1.1 of [23] and our Theorem 11.21 the expected number of sites 
visited by the walk in the first n steps and the expected area of the Wiener sausage over 
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the internal [0,n] coincide up to the error of 0(1) for the processes conditioned to return 
to the origin at the time n, provided that r is chosen to be r^, the variances of Brownian 
and random walk processes are the same and the fourth moments of the random walks 
exist. 

The results obtained in this paper are quite parallel with those in [2l], but the methods 
and the structures of the proofs are considerably different. In the random walk case we 
have simple expressions of the expected 'area' of the range of the walk for both free and 
conditional ones due to the discrete nature of the walk, which is not available to the 
Brownian case. In (cf. [2^) the Fourier analysis is the main tool, while in the present 
paper it plays a minor, though fundamental, role. We need a careful evaluation of certain 
integrals to find out the form of the second term as given in the formula of Theorem 11.21 
Our derivation of it is directed by the result of [21] , the agreement of asymptotic forms of 
various quantities for Brownian motion and Random walks being expected. The precise 
estimates of the first hitting time distributions for corresponding processes as obtained in 
[25] and [26] are used in both papers but the usage is much more essential for the present 
than in [25]. One can obtain an asymptotic estimate of the density of the first hitting 
time distribution valid uniformly with respect to starting positions and, with it, extend 
the result of Theorem II. 2l to the case when |x|/-\/t — > oo, which will be studied in another 
paper [27]. 

The following notation will be used: a/\h = min{a, b}, aVb = max{a, b} (a, 6 G R); two 
dimensional points are denoted by bold face letters z, x, y, |z| denotes the Euclidean length 
of z; and x ■ z the Euclidian inner product of x and z; for functions g and G of a variable 
X, g{x) = 0{G{x)) means that there exists a constant C such that \g{x)\ < C\G{x)\ 
whenever x ranges over a specified set; the letters C, C, C" etc. denotes constants whose 
values are not significant and may change in different places where they occur. 

We prove Theorem 11.11 in Section 2 and Theorem 11.21 in Section 4. In Section 3 we 
give some results on the distribution of the first hitting time to a disc, which prepare for 
the analysis made in Section 4. In the last section we derive the asymptotic expansions 
associated with N{t) as mentioned in Remark 1. 



2 Proof of Theorem 1.1 



We shall consider the Brownian motion started at z and denotes its law by Pz. Let 
a^'^^ = (Tu(r) be the first hitting time of Brownian motion Bt to the disc U{r) of radius r 
and centered at the origin. Then 



E[Area(Sf ^)] = / P,[a^'^ < t]dz. 



Let q{z;r) denote the density of the distribution of a^"^^: 

q{z,t-r) = jP^[a^^^ <t] 



so that 

d 



-E[km^{St-^)]= I q{z,t;r)dz. (2.1) 

at Jz>r 

Let pt{z) denote the heat kernel on the plane.: 



A 



As is well-known we have 



/ pt{z)e-^'dt = -Ko{\z\V2X) 

Jo TT 



(cf. HI, 0), hence 



/ q{z,t;r)e~^^dt 
Jo 



GaIOJzI) _ Ko{\z\V2X) 
Gx{0,r) Ko{rV2X) 



(2.2) 



where Gx denotes the resolvent kernel for the 2- dimensional Bessel process and Kjy the 
usual modified Bessel function of second kind of order u. Put 



From fl2.ip and 



m(t; r) = 

we have 

m{t] r)e~^^dt 



\z\>r 



q{z, t] r)dz. 



2nr^^^^udu 



27irKi{rV2X) 



KoirV2X)V2X 



(2.3) 



where we have applied the identity {d/dz)[zKi{z)] = —zKq{z) for the second equality. 
From the scaling property of Brownian motion it follows that 

m{t,r) = m{t/r'^] 1) 

and we suppose r = 1 in what follows. By Laplace inversion 



m{t, 1) 



-oo Kn{\J2iu)y2iu 



e'^'du. 



(2.4) 



Throughout the paper Igz and ^/z denote the principal branches in — vr < arg^; < tt of 
the logarithm and the square root, respectively. 

Although one can derive a leading term of m(t; 1) directly from (12. 4p as in [26] (the 
proofs of Lemmas 4 and 5), here we use the formula 



iV(t)e-A*rft = 1 _ 1 (A > 0) 

A — 1 A Ig A 

(cf. [11], page 196; also [7]), which somewhat simplifies the proof. Put 



(2.5) 



Then (12.30 is written as 



i^i(V2i) 



1 



/s:o(v2z)V2z K 



K ^Z —\ K ^Z Ifff/t ^z\ 



oo ^oo 

m(t; \)e-^^dt = 2n / N{Kt)e~^^dt + 2TT(p{X). 
Jo 



Accordingly (12. 4p becomes 

m{t, 1) - 27rN{Kt) 



(f{iu)e'^'^du. 



(2.6) 



Here (and in below) the trigonometric integral at infinity is improper. Note that f{z) is 
analytic on the slit domain — tt < arg^; < vr since Ko{^/z) has no zeros in it: the apparent 
singularity at z = k (i.e., that at A = 1 on the right-hand side of (12.51) ) is removable; 
also ip{iu) tends to zero as m — )■ oo and is bounded about the origin as will be observed 
shortly. 

The idea of the proof of the first formula in Theorem 11.11 would now be obvious. The 
asymptotic behavior of the Fourier integral on the right side of (12. 6 p for large values of t 
depends on that of (p{iu) near zero, provided that (p{iu) behaves sufficiently regularly. 

In view of the asymptotic formula 

1/2 / ^ \ 

7^,(^2^) = (^) e"^(l + ^^ + 0(l^r')J a« \^\^^ (2-7) 
(-TT < argz < IT, u >0) (cf. [H] (5.11.9)), we have 

^{iu) = ^=-^ + --^-— + R{u), (2.8) 

where R{u) = 0(|m|^^/^) with its derivatives R^^\u) = 0(|m|^^/^^'') as |m| — )■ oo. 
For the first formula of Theorem 11.11 it suffices to show that as t — )■ oo 



^ ^{iu)w{u)d'^du = ^(1^(1 + 0(1)), (2.9) 

where w{u) is a smooth function that equals 1 in a neighborhood of the origin and vanishes 
outside a finite interval, for the l — w{u) part contributes to the integral at most 0(l/t^) 
for any > 1 so that it is negligible. 
Put 

g{z) = - Ig (iv^ ) - 7 = -2-^ HK~h). (2.10) 

By definition 

MA-t^{ti-.-m.)) (2,11) 



and 



/ X 1 (^/2)"' z 1 ^ 1 



for —71 < axgz < n. It follows that for —1/2 < u < 1/2 

Ko{V2ki) = g{iu) + 2-Hu{g{iu) + 1) + Riiu); (2.12) 

Ki(y/2iu) 1 1 / , , 1\ „ , , , , 

1 1 —iu ( 1 \ , , , , 

'l + — ^ +i?3(M). (2.14) 



and 



Ko{\/2iu) g{iu) 2g{iu)\ g{iu 

Herei?i(M) = O^u^) x g(iu), R2{u) = 0{u)xg{iu) andi?3(M) = 0{v?) / g{iu). Substitution 
of these together with an easy computation yields 



where R^iu) = 0{u) with the derivatives R'^iu) = 0(1), = 0{l/ung\u\) {j = 

1,2). For evaluation of the contribution of R4 to the Fourier integral in (12. 9 p we split its 
range at m = ±l/t. The inner part plainly gives the bound 0(l/t^), whereas for the outer 
part J\u\>i/t RA{u)w{u)e^*'^du we repeat the integration by parts four times, which leads to 
the same bound of 0(l/t^) (see Lemma 2.2 of [25] for more details). Hence 

We rewrite the integral on the right-hand side as 



^ rioo ^K,tz 



i J-ioo [\gz] 



dz 



and apply the Cauchy integral theorem. The integral along the lower (resp. upper) half 
of the imaginary axis equals the one along the lower (resp. upper) side of the negative 
real axis in the positive (resp. negative) direction. Noting lg(— x ± iO) = Igx ± in for 
X > we then find the foregoing integral equal to 



00 



, -47r(lgx)e ^^"^ , Att . ... 

(as t — )• 00). Thus the relation (12. 9p . and hence the first formula of Theorem II. H has 
been proved. 

Consider the case t J, 0. Here we go back to the original inversion integral in (12. 4p . 
From the asymptotic formula (12. 7p we observe as before that the integrand of it involves 
the singular components 1/ \/2iu and 1/Aiu that are not integrable at infinity and we 
evaluate the contributions of them separately. To this end we bring in a function il){z) by 



IS ^ K^mz) 



Ko{^/2z)y/2^ v/2i 4(z + l)' 

so that i){z) = 0(1 z\ ^/^^ as 2; — 7- 00 in the domain — tt < arg z < tx where il){z) is regular. 
In addition we have il){z) = 0(1/2;) if 2; — and apply the Cauchy integral theorem to 
see that the principal value integral p.v. ip{iu)du vanishes, so that as t J, 



p.v. 



ij{iu)e''''du = / ^(w) e**" - 1 du = 0{Vi), (2.15) 



where the last equality may be verified by splitting the integral at m = ±l/t. 
On the other hand 

1 itu ' f°° COS u + sin u , 1271 



e'*''du= j= du=\—, (2.16) 



3 V2m ^0 ^/u V t 

and 

1 itu J 1 ["^ COStU + USmtU TX _^ 

— -e du = - / du = -e . (2.17) 

A{iu + 1) 2 Jo u^ + 1 2 ^ ^ 

From (I2.15P , fl2.16p and (I2.17P we find the formula as t J, as asserted in Theorem 11.11 
The proof of Theorem II. II is complete. 



/ 

J— ( 
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3 Preliminary results on the first hitting time to a 
disc 



Here we give several results which prepare for estimation of the expected area of the 
sausages for the Brownian bridges. Recall that g(z, t; r) denotes the density of the distri- 
bution of cr^^^ the first hitting time at the disc U (r) of Brownian motion Bt started at z. 
The following result is proved in [26] . 



Theorem 3.1 Uniformly for |z| > r, as t —> oo 

r 27lg(t/z2) , ^/ 1 



lg(|/t(z/r)^) _^2/2t ^ I 



/IN, i 0( /or z^ < t, 



(3.1) 



Let (T = cr^^) and q{z,t) = q{z,t; 1). We shall use the estimate of this theorem in the 
following slightly reduced form. 

Corollary 3.1 Uniformly for |z| > 1, as t — t- oo 

.(M) - W^2.M^)^0C-MM^) |zr<4Hg(lgt), 



OIttAtt) /or |z|^>4tlg(lgt). 



Proof. Using the inequality lg(z^/t) < z^/t for > t the assertion is immediate from 
Theorem 13.11 □ 

The following crude bound is often useful for dealing with the case |z| > i/GtlglgT. 
Lemma 3.1 There is a constant c > such that for all\z\> 1 and t > 1, 

g(z,t)<cpi+i(z). (3.2) 
Proof. For any unit vector ^, Ps{C) depends only on s and we see 

Pt(z) = [ q{z,t- s)ps{C)ds. 
Jo 

Hence 

Pt+i(z) > / g(z,t+l-s)— ds 

Jo ZTTS 

and, applying the parabolic Harnack inequality (cf. [20]) which implies that q{z,t) < 
Cq{z, t + 1 — s) for < s < 1, |z| > 2, we obtain the inequality of the lemma. (Note that 
for |z| < 2, a better estimate is given in Corollary 13.11 ) □ 

Remark 3. On the right-hand side of (13. 2p we can replace pt+i{z) by pt{z) for 
|z| < Mt, provided that at the same time c is replaced by a constant cm which may 
depend on M. It is warned that the constant cm actually depends on M and in fact the 
ratio q{z,t) /ptiz) tends to infinity whenever \z\/t — )• oo, t — )■ oo. 
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Lemma 3.2 For all |z| > 1 and t > 0, 

P^[a <t]<J^{ A l\e~^\^\-'^''^'. (3.3) 

V TT Y |z| — 1 J 

(2) 

Proof. The relation follows from the one dimensional result that if Bl denotes the 
vertical component of Bt, P(Q^\^\)[Bf'^ > 1 for < s < t] = 2(27rt)-i/2 j^oo^^^^ e-'^^/^^du. □ 



Lemma 3.3 There exists Ci > such that for |z| > 2 and t>l, 

d 

< J^E^\Ba ■ z; a < t] < CiyPt+i(z). 

Proof. The expectation EjJBfj ■ z; a < t] is a radial function of z and we may suppose z is 
on the upper vertical axis. Let tq be the first exit time from H \ ?7(1), where H denotes 
the upper half plane. Let bI^^ and Bf^ be the horizontal and vertical components, 
respectively, of Bt. Then by symmetry, for z = (0, ?/), y > 1, 

f^E^o,y)[B. ■z;a<t]= y^E(o,,) [4?; ^n, G H, tq < t]. 

The derivative on the right-hand side is expressed as the integral of the vertical component 
of ^ G H n dU{l) by P(o,y)[-BTo ^ d^^TQ G dt]/dt. It therefore suffices to show that 

|P(o,,)[i?.oeH,ro<t]<c|e-^^/^*. (3.4) 

For verification let D = H — (0,1) = {{x,y — 1) : {x,y) G H}. Then, on denoting by 
T£, = t{D) the first exit time from D, 

Piod^D e dt, -l< B^^^^^ < 1] 



> 



[ f ds-^P^o,y)[Bro e d^, To < s]P^[BI'^^. G (-1, 1), + S E dt]. 



Since ^ 

P{x,y)[rD e dt, edu]=pt{{x-u,y + l))^——dtdu, 

we obtain 

+ ! ^ /• ro / ,2pi_,(e+(0,l)) 



2pt((0,i/ + l))^ > / / -P(o,,)[5.„ Grfe,ro<s] 

t j|5|=i,CeHJo as t — s 

rt-l/A d 

> c — P(o,j,)[Pro G H,ro < s]ds. 

<^ t~~- 1 /2 Cio 



(is 



for some universal constant c > 0. Thus the bound of (13.41) follows in view of the parabolic 
Harnack inequality as before. 



n 



4 Wiener sausages for Brownian bridges 

For X G R^, put 

F{t, x;r)= f dz f f P^[(t, e ds, B,^ G rf^]pt-s(x - z - 0- 

J\x\>T Jo J\f\=r 



Then 



'|z|>r Jo J|5| = 

E[Area(5f'^) I 5^ = x] = ^t^(^> ^) + (4-1) 
From the scahng property of Brownian motion it follows that 

F(t,x;r) = F(t/r^x/r; 1) 
and we have only to consider the case r = 1 as in the proof of Theorem 11.11 Put 

Fo{t,x.) = [ dz [ Pj[a & ds\pt-si7--^) = 1 dz j q{z, s)pt-siz - x.)ds. 

J\z\>l Jo J\z\>l Jo 

We first consider the case x = in Subsection 4.1. The general case, dealt with in 
Subsection 4.2, is based on the proof of this special case but need additional estimations 
of a certain integral that are involved and partly delicate. 



4.1 The case x = 

Lemma 4.1 Let F and Fq be as above. Then F{t,0; 1) = Fo(t,0) + 0(l/t). 

Proof. Let a be a constant larger than 1. In this proof it may be arbitrarily fixed (eg. 
a = 2; in the case x 7^ treated later we shall take a = x^). We split the range of 



integration by the surfaces s = t — a and |z| = -y/4(t — s) \g{t — s) put 
Dq = DQ{a) = {(s, z) : t - a < s < t, |z| > 1} 



D> = D>{a) = {{s,z) : < s <t - a,\z\ > JA{t - s) \g{t - s), |z| > 1} 



L)< = D<(a) = {(s, z) : < s < t - a, |z| < ^ 4:{t - s) \g{t - s), |z| > 1}. 

Accordingly we break the rest of the proof into three parts. Some estimates obtained 
below are more accurate than necessary for the poof of the present lemma, but they are 
needed in the proof of Lemma 14.41 essential for our proof of Theorem II. 2[ 

Part 1 -.Dq. The contribution to the integrals defining F and Fq from Dq is 0(l/t(lgt)^). 
Indeed on the one hand we can use the bound (13.21) for the integration w.r.t. z over the 
range |z| > yjAt Ig Igt, in which pt+ii^) = 0{l/t(\gt)^). On the other hand, on applying 
Corollary 13.11 the contribution of the other part is at most a constant multiple of 

/ -7^ ^ / (lg|z|) SUpPf_,(z-Oc?Z 

Jt-a S{igSy Jl<\z\<y/mglgt 1^1=1 



1 n 



Part 2 :D>. For (s,z) G D^, according to Corollary 13.11 q(z. s) < C/s(lgs)^ if 
a < s < t/2 and g(z, s) < C/s\gs if s > t/2. Hence both the contributions of D> to F 
and Fq are dominated from above by a constant multiple of 



t—a 



l(s < t/2) l(s > t/2) 



< 



t2 

o 



S(lgs)3 

c 



C rt-a 



s Ig s 

ds 



ds 



|z|>^4(t-s)lg(t-s) 151=1 



sup Pi-s(z - ^)(iz 



t/2 (slgs)(t-s) 



t\gt 



On using fl3.3p the integral on < s < 4 is readily evaluated to be at most 0{l/t^). 
Part 3 :D<. Observe that if (s, z) G -D<, then 

O-pt-.(z) = pt_.(z)(e[^-^-^l/(*-^)-l) 



with 



Vt{s,z,^) =pi_^(z) X O 



t 



+ Vt{s,z,^) 



(4.3) 



z^ + l 



The integral for the difference F{t, 0; 1) — FQ{t, 0) restricted to 

Z^o := {(s,z) G D< : < s < 4} 

is at most 0(1/^^) in view of ( 13. 3p . Thus we may restrict the integral to the range D^\D^. 

Consider the contribution of the first term on the right side of (14. 3p . Applying Lemma 
13.31 with the help of Corollary 13.11 (for |z| < 2) we deduce that 



dsdz 

D<\D% 

rt—a 

<C ds 

J4 J\x\>2 



z-^ P^[a G ds,B„ G d^] 

Pt-siz)- 

151=1 t-s 



ds 



Pt-^s{z)z^Ps+l{z) 

{t - s)s 



ds 



4 (t-s)2s(lgs)2 - t 



< ^- (4.4) 



Here we have applied the trivial bound pt-si^z) < 1/it — s) (for the integral on s < t/2) as 
well as Ps+i(z) < 1/s (for that on s > t/2). The part involving the term 1/2 is evaluated 
to be 0{l/t) by dominating q{z,s) by Cps+i{z) for s < t/2 and by C/slgs for s > t/2. 
For the proof of the lemma it now suffices to show that 



Jd^\d'^ 4|=i ' ^ ds \t\gt 



(4.5) 



To this end we further split the region D< by the surfaces s = t/2 and |z| = 1 + 
y/As Ig Ig s (4 < s < t/2). Denote by J^"(a), , the double integrals on the 
regions 



D 



RH, 



< 

:= {is,z) eD^:A<s< t/2, |z| > 1 + ^4slglgs }, 
:= {(s, z) G D< : 4 < s < t/2, |z| < 1 + v'4slglgs }, 



a 



{(s,z) G D< : t/2 < s < t - a}, 



respectively. Employing Corollary 13.11 we obtain 

rt—a rj a 

J^^fa) < C ^ 



t/2 slgS y|z|<^4(i-s)lg(t-s) it-s)' 
t Ig t J t/2 t — S \t J 

and, similarly but dominating pt-si^) simply by — s), 



4S' < c 



4 {t- S)3 y|z|<l+^4slglgs Ig S 

< -7 / — ds = 0[ — — 
Ja Igs Vtlgt 



Ps{z)dz 



With the help of ([32]) and (O) as well as g^]) we deduce that 



t/2 



./|z|>l+^6slglg; 
+ 



z^Pt_, Z 
Ps+i(z)^^ —dz 



it - sf 

Psiz) (Iglgs)^] z2pi_,(z) 



+ 



|z|>l+^4.1glg. I IgS Z2(lgs)3j {t-S 



dz 



.t/2 

< C ds 



+ 



z'^Ps{z)dz 



+ 



{\gs)H^ J\z\>l + y/4slglgS t^lgS J\z\>l + y/6slg\gS 



< 



C" 



Thus (14. 5p has been proved. The proof of the lemma is complete. 



Lemma 4.2 



Fo(t, 0) = NiKt) - nptil) + 0(l/t(lgt)2). 
Proof. The Laplace transform of Fq is given by 

/ Fo{t,0)e~^*dt = / e'^^dt dz P^[a E ds]pt^siz) 

Jo Jo J\z\>l Jo 



Koi\z\V2X) 



TC J\z\>l 

(cf. [12] Sect. 7. 2). Fortunately we have the identity 

1 d r 



dz ■ 



(Klira) - Kl{ra)) 



for any constant a > 0, so that 



^0 Kq 



[i^i(V2A)] 



Hence 



Fo(t,0) = -7rpt(l) + — 



277 J-oo ii'o(v2m) 



2A) 



e'^^'du. 



For evaluation of the last integral we can proceed as in the preceding section. Put 



Eiz) 



[K^iV2-z)l _ 1 



1 



1 



K — 1 K ^Z lg(K ^z) 



Then for -1/2 <u< 1/2, E{iu) = -1 + i + [Ag{iu)]-'^ - [2g{iu)]-^ + R{u) with R{u) 
0{u) and as before (see (12.51) and the ensuing discussion up to (12. 6p ) we have 



Foit,0)-N{Kt) + npt{l] 



The proof of the lemma is complete. 



2tt 



E(iu)e''''du = O 



□ 



Combining Lemmas 14.11 and 14.21 with (14. ip leads to the assertion of Theorem 11.21 with 
x = 0. 



4.2 The case x 7^ 0. 

We first extend Lemma [4. 1[ 

Lemma 4.3 For each M > 1 uniformly for |x| < My/t, as t —t- 00 

F(t,x;l) = Fo(t,x)+0(l/t). 
Proof. First of all we verify that 

/ rfz /* P,[a e ds, G rfeh-s(z - X - < 7^. (4.8) 

J|x-z|<l,|z|>l Jt-l t\gt 

For verification we use the inequality 

Pz[o- G ds, edE\< j dypi{y - z)Py[a + 1 E ds,B„ E d^]. 

If |x — z| < 1, pi(y — z) is dominated by e^^'^^^^^''^. Hence, performing the integration by z 
first, we see that the repeated integral in (14.81) is at most //_^ ds J e~^'y^^^'^^'^q{y, s — l)dy, 
which is evaluated to be 0{l/tlgt) according to Corollary 13.11 

The rest of the proof proceeds in parallel with that of Lemma 14. 1[ Of course we must 
replace pt-s(z — ^) and pt_s(z) by pt_s(z — x — ^) and by pt-si,'^ — x), respectively. We 
accordingly modify the definition of D< as 



L)< = {(s, z) : < s < t - a, |z - x| < ^4(t - s) \g{t - s)] 

and analogously. Then Parts 1 and 2 are treated without any change except that for 
Part 1 we use (14. 8p . For Part 3 we consider 

Pt„,(z-x-0-Pt-.(z-x) = j9t_,(z-x)(e[(^-^)-«-^]/(*-)_l) 

= p,_,(z-x)^^^i^^-i^+r/i,,(s,z,0 (4.9) 

b S 

with ^ 

^t,x(s,Z,0 =Pt_.(z-x) X ^J^2 ^ 



in place of (14.31) . In the first term on the last member of fl4.9p x ■ ^ may be replaced by 
(x ■ z)(z ■ ^)/z^ since by symmetry the component of ^ perpendicular to z vanishes after 
integration. It follows that (z — x) ■ ^ can be replaced by (z ■ ,^)z ■ (z — x)/|zp. Keeping 
this in mind and employing (13. 3 p together with Lemma 13.31 and Corollary I3.lt we deduce 
that 



/ dsdz 


/ Pt- 






C 1 


ft— a r 




ds / 




2 J\2 




f dz 1 




|z|<2 J 



Pt-s{Z - Xj- 

|z ■ (z 



- x) ■ ^ - I Pr,[a e ds,B^ e d^] 



it 



S]S 



X e 



- s ds 
Pt-siz - x)psiz)dz 

^|z-x|2/2(t-s) 



{t-sys{\gs) 



< 



c 



(4.10) 



Here the last inequality is due to the observation that the first integral in the middle 
member restricted on the interval s G [a,t/2] is at most a constant multiple of 



1 /•*/2 , /• z^ + |z||x| , , , 1 Ixl rVa i /i 
— / ds ps[z)dz<- + — —ds = 0[ — 

J a JR2 S t t'^ Jo y/s \t 

and similarly for the other interval. 

The rest of the proof is similarly dealt with and hence omitted. 



□ 



Lemma 4.4 Foe each M > 1 uniformly for |x| < M\/t, as t ^ oo 



F„{t,xj - F„{t,0) = 2ir[p,(x)-p,(0) 



t X^ 

+ 



Igt 2t(lgt)^ 



t 



X2 V 1 



+ 0(1) 



tlgt 



The right-hand side of the formula of this lemma is suggested by that of Theorem 11.21 
and the latter in turn is by the corresponding one in [2l] for random walks. The course 
of proof given below is often steered by the asserted formula. 

Proof of Lemma \4^\ We suppose |x| < \/i /2 for simplicity of the description of the 
proof (otherwise one may take a = x^/2M^ instead of a = x^ in below). Suppose also 
|x| > 2, the same proof of Lemma 14.11 being applicable with little alteration if |x| < 2 
. We make computations analogous to those in the proof of Lemma 14.11 We substitute 
a = x^ in the definitions of Dq, Dy and D< given therein and denote the resulting regions 
by Do(x^); -D>(x^) and D<(x^), respectively; denote the contributions from them to the 
difference Fo(t,x) - Fo(t,0) by /[DqIx^)], /[£'>(x2)] and I[D^{x.% respectively. 

For the estimate given in Part 2 of the proof of Lemma 14.11 where the integral over 
D>(a) is dealt with we have the same bound 0(l/tlgt), namely 



j[D>(x2)] = o(iAigt), 



(4.11) 



since |z — x| > |z| — |x| > v^3tlgt if t/2 < s < t — x^ and (s,z) G D>(x^) (the integral on 
< s < t/2 is dealt with in the same way as before). 



The rest of the proof is somewhat involved and broken into three parts 1 to 3. 
1. Estimation of I[Dq{'x.'^)]. There exists a constant C such that for |x| > 1, 

Pu{z)dz < Cx^ 







r du [ 


Ig — Ig z — x^ 


Jo JR2 





(4.12) 



as is proved shortly. For the integral that gives /[i5o(x^)] we may restrict (s,z) to {4 < 
t—s < x^, |z| < a/4s Ig Igs}, the contribution of the remainder being at most 0(x^/t(lgt)^) 
as discussed before. In view of this as well as of f l4.12p we may write it as 



t-4: 



ds 



x2 y|z|<^45igig3 (Igs)^ 
+ 0i 



(Ig |z - xHpt_,(z - x) - (lgz2)pt_^(z) 



dz 



Substituting from the identities Ig |z — xp = lg(t — s) + lg[|z — xp/(t — s)] and Igz^ = 
lg(t — s) + lg[z^/(t — s)] and noting that /lg(z^/(t — s))pt-siz)dz equals /(lgz^)pi(z)(iz, 
a finite constant we find that 



/[Do(x')] = 27r pi(x)-pt(0) 



-4 lg(t - s)ds 

t-x2 (Igs) 2 



+ 



t{\gty 



(4.13) 



This is not negligible and to contribute to the leading term on the right-hand side of the 
formula of the lemma. 

Proof of fl4.12p . If the range of integration is restricted to |z| A |z — x| > |x|/2, then 
Igz^ — lg|z — xp is uniformly bounded and the corresponding (repeated) integral is 
obviously bounded by a positive multiple of x^. The integral corresponding to {|x — z| < 
|x|/2} can be readily evaluated to be ©(x^). For the range |z| < |x|/2 it suffices to show 
that 

X^ Jo U J|z|<|x| V \z\ ) 

Changing the variables of integration according to s = ujy? and y = z/|x| transforms 
the left-hand side of (14.141) to f^s'^ds /|y|^]^(lg |y|)e"^ ^"^^dy, which is independent of x. 
By another change of variables this last integral is further transformed to 



°° ds 



/I s^ J\y\<^s \ \y\ 

which is certainly finite. Hence we have fl4.14p . 

2. Estimation o/J[D<(x2)]. For (s, z) G D^i^), instead of (14. 3 p we have 



Pt~s{z - X) -pt-s{z) 



Pt_,(z)(e[^-^-^^']/(*-^)-l) 

(z ■ x)2 (z ■ x)x2 



1 2 

Z ■ X — pC 



with 



t - s 
+ ?7f,x(s,z) 



?7t,x(s, z) = pt-siz)-^^ X O 



+ 



2{t-sy 2{t-sy 



(4.15) 



t-s 



t-s {t-sy 



Since the double integral of the terms that are linear in z inside the big square brackets 
in fl4.15p vanishes by skew symmetry, the substantial part of pt_s(z — x) —pt-s{z) reduces 
to pt_5(z) X Li_x(z), where 



^t,x(s,z) 



+ 7rr. -T^ + - X 0{- + 



2(t - s) 2{t - s)2 t-s 



t- s (t- s) 



(4.16) 



Let DO , L'™(x2), and D^^ be defined as in Part 3 of the proof of Lemma 13.31 
and denote by /[-D<], /[-D^^(x^)] etc. be corresponding double integrals on them, so that 

/[D<(x^)] = + /[^^"(x^)] + I[Di^] + /[D^]. (4.17) 



The integrals I[D^], I[D^] and /[Dlh-'] are independent of x as the notation suggests. 



2.1. Estimation o//[D^^(x^)]. By definition 



/[/^^^(x')] 



ds / g(z,s) 

t/2 J\z\<y/4{t-s)lg{t-s) 



dz. 



(4.18) 



Let (s, z) be in the region of integration of this integral. By what is remarked right 
above, pt-s(z — x) — pt-s{z) may be replaced by pt_s(z)Lt^x(s, z), which is bounded in 
absolute value by a constant multiple of 



Pt^siz){t- S)-'^' 1 + it-S 



Keeping this in mind one replaces the upper limit i/4(t — s) lg(t — s) of the inner integral 



in (14181) by yJ8(t - s) lglg(t - s). Since g(z, s) < C(lg |z|)/s(lgs)^ in view of CorollaryO 

and the integral of (Ig |z|)[l + (t — s)^^z^]]9j_s(z) w.r.t. z over |z| > ^Js(t — s) lglg(t — s) 
is at most l/[lg(t — s)]^, the error given rise to by this replacement is bounded above by 

rt-^^ (7s / x^ 



It/2 S{\gsm-S){\g{t-S)y 

This together with Corollary 13.11 shows that 



/[^^"(x^)] 



ds 



Igz 



t/2 J\z\<y/8(t-s)lglg{t-s) (lgs)2 

2 

o' 



2vrp,(z) 



t(lgt)^ 



Pt-s[Z - Xj -pt_,(zj 



t(lgt)2. 

Now, substituting the decomposition Igz^ = lg(t — s) + lg[z^/(t — s)], we have 

*-x2 lg(t - s] 



I[Df'{^')] = 2n[p,{^)-p,{0) 
where 

ir(x,t) 



t/2 {\gsy 



-ds + K{x,t) + 



tilgty 



dz 

(4.19) 
(4.20) 



t-^' ds 



t - s 



t/2 (lgs)2 y|z|<V8(t-s)lglg(t-s)) 

On writing m for t — s, using the identity 

2TTPu{z)pt^u{z) = t~^P(t-u)u/tiz), 



27rps(z)pt_s(z)Lt,x(s, z)dz. 



noting that one may replace \g{7? /u) by Ig {tP' /\u{t — u) /t]^ in the present estimation, and 

changing the variables of integration according to y = z/ ^JT(T-^'s)Jt, the inner integral 
above becomes 



lg(zV«) 



\z\<y/8u IglgM t 



P(t-u)u/t{2)Lt^x{t - u, z)dz 



igy 

\y\<y/8{t/{t-u))lg\gu t 



2 r 



x2 (t-M)(y-x) 



2n 



1 7 ^ ~ ^7 \ 2 ^ 



2m 2M)f: 

2 

^ +R 



pi{y)dy + R 



tu{\gufl^ 



(under the restriction < u < t/2), where R denotes the contribution of the error term 
in Lt^y(_{t — u,z) and we have used the identity /(y-x)^/(|y|)(iy = 2~^x^/ y^/(|y|)(iy valid 
for a non-negative function /(r) and 

bi= (igy^)pi(y)c?y and &2 = / y^(igy^)pi(y)c^y- 

It follows that 62-261 = 2and \R\ < C\K{x,t)\x^ /t. Noting that J^i'^ u-\\g{t-u))-'^du = 
(lgt/x2)(lgt)-2 + 0((lgt)-2) we infer 



K(x,t) 

and substitution into f l4.20p yields 



/[D^«(x2)] = 2vrL(x)-p,(0) 



x2 lg(t/x2 

2t(lgt)2 



+ 



2 / ' 



t/2 (Igs)^ 



2t(lgt)^ 



t(lgt)^ 



(4.21) 



2.2. Estimation of I[d\^{'x'^)]. This part is similar to the preceding one but much 
simpler. This time we substitute Igz^ = Igs + lg(z^/s) to see that 



/[4<H^(X^)] 



t/2 



ds 



|z|<^4s Iglgs 



q{z,s) 



Pt-s[z - xj -pt-s[z} 



dz 



27T 



Ptix)-ptiO)] 



*/2 ds 
4 IgS 



+ R 



with 



Ja le s r 



< 



M (Ig s)2 JR2 " 
C"x2 



Z2x2 



t- s (t-s) 



Z^X 

Ig — )Ps(z)c?Z 



(4.22) 



t(lgt)2- 

2.3. Estimation of I[Dl^{x^)] and /[D^h^x^)]. It is readily checked that (x^)] 



0(l/t(lgt)^). By the same argument as made for the estimation of J^J^ in the proof of 
Lemma [4. II we deduce that 



/[^S(x')]=0 



(4.23) 



1 T 



3. Completion of Proof . Combining KT^ . dOT]) . ( lO^ . (lO^ and ( HTT]) we deduce 



that 



J[D<(x2)] = 27rpi(x)-pi(0) 



lg(t/x^ 
' 2t(lgt)2 



+ 



t/2 (lgs)2 
,2 



*/2 rfs 

Igs 



t{\gty 



Finally, noting ^pfip^t^s + /o^^ = V Ig^ + 0(t/(lgt)^) and recalling fHTT]) as well as 
Fo(t, x) -Fo(t, 0) = I[Doi^^)]+I[Dy{x^)]+I[D<{x.^)] we conclude the formula of Lemma 

sai □ 



Proof of Theorem \l.^ From Lemmas 14.21 and 14.31 it follows that for 1 < |x| < M\/t, 
Fit, X, 1) N{Kt) - 7rpt(l) + Fo{t, x) - Fo{t, 0) 



By Lemma 14.41 we can write the right-hand side as 



+ 0(1). 



2TltN{Kt) + 2TXt 



N{Kt) + 



Igt 



Pt(x) -pf(0) vrx 



2 r 



+ 



(Igt)^ 



1? 



+ 0(1) 



+ 0(1), 



in which the second term reduces to 0(x^/(lgt)^) since N{nt) = 1/ Igt + 0(l/(lgt)^). In 
view of (14. ip this yields the formula of Theorem II. 2[ □ 



5 On the asymptotic expansion related to A^(A) 

Our arguments made in this section are based on the following formula due to C. J. 
Bouwkamp [2]: for A > 0, a > 

du= / r(cr + x)e ''^^^dx+ ^ ' \ (5.1) 



ilguY + n"^ Jo vr vr^A 

where s may be an arbitrary positive number and |^(A,s)| < 1 (Eq (11) of [2j; see also 
Doming et al [1] Eq (17) for the case a = 0). From this he obtains the asymptotic 
expansion of A^(A) in powers of 1/ Ig A as A — ?■ 00. We shall extend the argument of [2]. 
Take o" = in (15. ip so that the integral on the left gives A^(A) and substitute from Euler's 
relation r(x)r(l — x) = n/ sinvrx as well as A = at to see that as t — ?■ cxd 

N{at)= V e-^'g^rfx + of-V (5.2) 

^ ' Jo r(l-x)a^' VtV ^ ' 

Lemma 5.1 For any constant a > 0, the function N{at) admits the asymptotic expan- 
sion 

N{at) ~ E [t ^ 00) 

^ ^ (lgi)„=o(lg^)" 

with the constants ««(«) determined by 

1 (1^1 <1). (5.3) 



r(l - x)a^ ^0 n! 



Proof. The integral on the right-hand side of (15. 2 p is a Laplace transform, in the variable 
Igt, of a function regular at a; = and it is standard to infer that the integral admits an 
asymptotic expansion in the powers of 1/lgt with the coefficients determined by (I5.3p . 
showing the lemma. □ 



Lemma 5.2 For each constant a > 0, as t ^ oo 



N{as)ds 



t 



E 



bJa) 



(lgt)„t^,(lg^)"' 



— asN'(as) 



ds 



E 



(n + l)6„(a) 



with the constants 6„(q;) determined by 

1 



oo t 

E- 



{l-x)T{l-x)a' 
Proof. Integrate the both sides of (15. 2 p to see that 



< 1). 



l/a 



N{au)du 



t 



-X let 



x)T{l — x)a 
l/a 



dx 



'l-x)T(l-x] 



dx + O 



(5.4) 



which is valid for every s > 0. Note that 1/(1 — z)r(l — 2;) = l/r{—z) is an entire function. 
Since 0(1) term is negligible for the asymptotic expansion, we obtain the first formula as 
in the preceding proof. One can verify the second one in a similar way but by employing 
(15. ip with (7 = 1. Alternatively, note that the left-hand side integral of the second formula 
equals Jq N{as)ds — tN{at). Then, we readily derive it from the first one combined with 
the preceding lemma. □ 



From the Weierstrass product formula 

1 



r(i + x) 

(cf. [H]) it follows that if ({z) = J:T=i 

1 



h 



r(i -x) 



00 2 

-7a; - E -C(^^)a;", 



n=2 



n 



hence 



E 

n=0 



ni 



-X = exp 



00 

(7 + lga)x- ^ -C(^)a;" 



n=2 



n 



Similarly, using (1 — x) ^ = exp J2'^=i n ^a^" "we obtain from (15. 5p 

h (fy) r °° 1 

E = exp - (7 + lg« - l)x - E -iCH - l]x^ 



n=0 



n=2 



n 



(5.5) 



From these identities one can readily computes the first several terms of (a„(a;))^Q and 
those of (6„(a))^g as exhibited in (II. 2p and (II. 3p . respectively. 



1 n 
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